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Collapse calculations indicate that the hot newly born protoneutron stars (PNS) rotate differentially so that strong toroidal 
magnetic field components should exist in the outer crust where also the Hall effect appears to be important when the Hall 

parameter f3 = lobt is of order unity. The amplitudes of the induced toroidal magnetic fields are limited from above by 
the Tayler instability. An important characteristic of the Hall effect is its distinct dependence on the sign of the magnetic 
field. We find for fast rotation that positive (negative) Hall parameters essentially reduce (increase) the stability domain. It 
is thus concluded that the toroidal field belts in PNS induced by their differential rotation should have different amplitudes 
in both hemispheres which later are frozen in. Due to the effect of magnetic suppression of the heat conductivity also the 
brightness of the two hemispheres should be different. As a possible example for our scenario the isolated neutron star 
RBS 1223 is considered which has been found to exhibit different X-ray brightness at both hemispheres. 
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Progenitors of neutron stars are high-mass stars with more 
than eight solar masses that develop a degenerate iron core. 
If the core mass approaches the Chandrasekhar limit it be- 
comes gravitationally unstable and implodes. The collapse 
comes to a temporary end if nuclear densities are reached. 
At that stage the rebounding inner core drives a shock wave 
into the outer core, a mechanism that is currently believed 
to be responsible for the appearance of supernova. 

If the core of the supergiant rotates already rapidly the 
neutron star will be born as a fast rotator with an angular 
velocity near the break-off value, i.e. 1 kHz. This value ex- 
ceeds the rotation rate of the fastest young pulsars known 
by one order of magnitude so that the question arises how a 
critically rotating protoneutron star (PNS) spins down. One 
possibility is the angular momentum loss by gravitational 
wave emission via unstable r-modes (Friedman & Schutz 
1978; Andersson 1998; Stergioulas & Font 2001; Lindblom, 
Tohline & Vallisneri 2001). As the viscous damping of the 
r-modes is smallest at temperatures around 10° K, this insta- 
bility works best as long as the neutron star remains hot. Up 
to 90% of the rotational energy can be removed in that way 
from the newly formed neutron star within hours. Other pos- 
sibilities involve angular momentum transport due to non- 
axisymmetric instabilities also connected with gravitational 
waves. 

Following Burrows (1987), entropy-driven convection 
may play an essential role in the neutrino-mediated super- 
nova explosion scenario since it enhances the neutrino lu- 
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minosities in the post-collapse stage. Such a convection 
might be important with regard to a standard-dynamo ac- 
tion in PNS (Thompson & Duncan 1993). If, on the other 
hand, additionally differential rotation exists in the turbulent 
domain, then an a!7-dynamo can work producing strong 
toroidal magnetic fields (Bonanno et al. 2003, 2006). In- 
deed, hydrodynamic simulations of rotational supernova 
core collapse have shown that even a nearly rigidly rotat- 
ing initial core results in a strongly differentially rotating 
post-collapse neutron star (Monchmeyer & Miiller 1989; 
Janka & Monchmeyer 1989; Dimmelmeier, Font & Miiller 
2002; Kotake et al. 2004; Ardeljan et al. 2005; Burrows et 
al. 2007). Any nonhomologous collapse creates necessarily 
some degree of differential rotation if angular momentum 
is conserved locally during collapse. The models reveal a 
strong differential rotation in the azimuthally averaged an- 
gular velocity (Ott et al. 2005). Differential rotation may 
furthermore be generated by r-modes via nonlinear effects 
(Rezzolla, Lamb & Shapiro 2000) or simply by accreting 
falling-back material (Watts & Andersson 2002). 

1.1 Differential rotation and magnetic fields 

In the presence of a poloidal field Br, differential rotation 
with a shear q produces a toroidal field component by in- 
duction. The ratio of the resulting field to the original one 
can simply be estimated as 

e=-^~qttr if e < Rm, (1) 
Br 

with Rm as the magnetic Reynolds number of the differen- 
tial rotation. For high Rm the differential rotation may in- 
duce strong toroidal fields. Also flux compression will play 
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an important role in amplification of both poloidal fields and 
toroidal fields (Burrows et al. 2007). However, the result- 
ing magnetic field transports angular momentum outwards 
and feedbacks the differential rotation. The timescale of this 
backreaction is 

Ho p qSIL 2 po 9 qfJL 2 



BrB^ 
With e ~ qflT one finds 
V9o9 L 



D 



(2) 



(3) 



R 



i.e. the Alfven travel time of 1 ... 10 s (Shapiro 2000). 
Hence, after e ~ 100 . . . 1000. Typical values of the 
neutron stars have been used: p ~ 10 13 g/cm 3 , Q ~ 100 
s _1 , L ~ 10 5 cm (the crust thickness) and B r ~ 10 12 G. 
Note that for e ~ 1000 the differential rotation is immedi- 
ately destroyed. We find e < 1000 (i.e. B toI < 10 15 G) as a 
necessary condition for the existence of differential rotation 
over several rotation periods. This value is a rather small 
value insofar as 

HL 2 10 12 cm 2 /s 



Rm = 



(4) 



T) Tj 

so that already for rj < 10° cm 2 /s the critical e is exceeded. 
The microscopic value is only r\ ~ 10~ 6 cm 2 /s. It is obvi- 
ous that with such a small microscopic value a differential 
rotation cannot survive. It is an open question whether such 
high values of r\ can be reached in neutron stars (see Naso 
et al. 2007). 

The diffusion time L 2 /rj for rj ~ 10 9 cm 2 /s is also 10 s 
which, however, would also be the decay time of the differ- 
ential rotation for magnetic Prandtl number Pm > 1. With 
such high values of viscosity a prescribed differential rota- 
tion cannot exist longer than a few rotations. 

In the present paper we assume that differential rotation 
exists for at least 10 s (~ 100 rotations). During this time 
r-modes are excited producing gravitational waves. The vis- 
cosity must thus be < 10 9 cm 2 /s. This value corresponds to 
the expression v ~ ah 2 Q with a ~ 10~ 3 which is known 
from the accretion theory (due to small-scale MRI). As the 
microscopic magnetic Prandtl number Pm is very large for 
neutron stars we also continue with a high magnetic Prandtl 
number for the unstable PNS crust matter, say Pm ~ 100 
so that rj ~ 10 7 cm 2 /s. In that case the magnetic diffusion 
time exceeds the viscous time by a factor of 100 and the 
magnetic Reynolds number is of order 10 5 . The differen- 
tial rotation would thus generate huge toroidal fields with 
e ~ 10 5 which, however, would destroy the differential ro- 
tation. 

Therefore, in the present paper the stability of strong 
toroidal magnetic fields against nonaxisymmetric perturba- 
tions is probed in order to find their real upper limits. We are 
thus considering the Tayler instability under the influence of 
differential rotation and for high magnetic Prandtl numbers. 
The toroidal field is assumed to dominate the poloidal field 
(e > 1) so that stability under only toroidal field is con- 
sidered. To that end, as will be shown in the next section, 
also the influence of the Hall effect in neutron stars must be 
taken into account. 



1.2 Magnetic fields and Hall effect in neutron stars 

Neutron stars have the strongest magnetic fields found in the 
Universe, with fields exceeding 10 13 G for young (~ 10 7 yr) 
radio and X-ray pulsars, and a still appreciable 10 8 -10 10 G 
for much older (~ 10 10 yr) millisecond pulsars. This cor- 
relation between field strength and age suggests that these 
very different strengths are due to the field decaying in time 
rather than to differences between different neutron stars. 

Jones (1988) and Goldreich & Reisenegger (1992) have 
proposed that correlation between the magnetic field of the 
neutron star and its age is due to the Hall drift. Since the Hall 
effect enters the evolution equation for B as a quadratic non- 
linearity, it necessarily leads to a timescale inversely pro- 
portional to |B|. The Hall effect is therefore attractive for 
explaining the variations in the decay rates (for B ~ 10 13 
G the field should evolve on a 10 7 year timescale while if 
B ~ 10 10 G it should evolve on a 10 10 year timescale). 

There is a bulk of literature about the existence of the 
Hall effect in neutron stars. The main findings may be sum- 
marized as follows. The Hall effect strongly depends on the 
magnetic field amplitude and the temperature of the neutron 
star. In the presence of strong magnetic fields the magnetic 
diffusivity is anisotropic and is given by a tensor whose 
components along the magnetic field are 7711 , the components 
perpendicular to the magnetic field are r)±, and off-diagonal 
Hall component rjs. For more details concerning the gen- 
eralized Ohm's law in multi-component plasma we refer to 
the papers Yakovlev & Shalybkov (1991) and Shalybkov & 
Urpin(1995). 

With Hall effect included the magnetic induction equa- 
tion takes the general form 

<9B 

~dt 

with r/ = rj± and (3 = c/47ren e , where n e is the electrons' 
number density. In addition, it is useful to define the Hall 
magnetic diffusivity as r/n = (3B. The Hall effect becomes 
important if (3 > 1, where 

m 
v±' 

For magnetic fields smaller than some critical value, B CI , 
Vi- — r l\\ — ^0 where r/o is the magnetic diffusivity without 
an applied magnetic field. If B > B CI than rj± increases 
as B 2 for increasing magnetic field. The Hall magnetic dif- 
fusivity, on the other hand, is proportional to the magnetic 
field value. As a result, the Hall effect can be important only 
in some vicinity of the B CI . 

The critical magnetic field can vary significantly within 
the neutron star envelopes depending on chemical composi- 
tion, temperature and density. According to Potekhin (1999) 
the critical magnetic field is ~ 10 12 G for iron composi- 
tion with temperature 10 8 K and density 10 11 g/cm 3 . De- 
tailed calculations of the electrical conductivity in pure neu- 
tron star crusts (Cumming et al. 2004) indicate that the Hall 
time scale under such parameters is indeed shorter than the 



r/AB = curl(u x B - (3 curlB x B) 



(5) 



(6) 
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Ohmic decay time, which means that f3 > 1. One finds3 
that for iron (Z=26, A=56) with p = 10 13 g/cm 3 the ft 
varies from lCT 3 .Bi2 for T = 10 10 K to 3Bu for 10 8 K. 
Note, however, that for the same plasma the Hall parameter 
(3 reaches a maximal value of ~ 10 for B ~ 10 13 G and 
decreases for higher magnetic field values. 

Hence, it makes sense to ask for the consequences of the 
Hall term for young neutron stars with fields of Byi ^S> 1 
which can be imagined - and this is the point here - as 
toroidal field due to the induction of a differential rota- 
tion. Important for us is only the assumption that the (early) 
phase of the existence of differential rotation in the crust 
of the PNS is accompanied by (3 of order unity for the re- 
sulting toroidal fields. Obviously, (3 linearly depends on the 
magnetic amplitude so that we can write 

(3 = f3 VPmHa = f3 S, (7) 

with S as the Lundquist number (see below). The parameter 
f3o does not depend on the magnetic fielc0. 

To estimate the magnetic Prandtl number we should also 
use r/± instead of rj . We will have smaller magnetic Prandtl 
numbers for the parameters where the Hall effect is impor- 
tant. Nevertheless, it is easy to estimate that the magnetic 
Prandtl number can be much larger than 1 for the typical 
neutron star envelopes parameters. 

2 Tayler instability 

Differential rotation leads to an increase of the toroidal com- 
ponent of the magnetic field by winding up the poloidal field 
lines. The ratio of both components is given by the magnetic 
Reynolds number of the differential rotation which is very 
large also for PNS due to their small values of the magnetic 
diffusivity rj. The same is true for the microscopic magnetic 
Prandtl number (v ~ 10 cm 2 /s, r\ ~ 10~ 6 cm 2 /s). 

Too strong toroidal fields, however, become unstable 
against the Tayler instability. Tayler (1961, 1973) and Van- 
dakurov (1972) considered the stability of nonaxisymmetric 
disturbances and showed that for an ideal fluid the necessary 
and sufficient condition for stability is 
d 



dR 



{RBi) < 0. 



(8) 



An almost uniform field would therefore be stable against 
axisymmetric perturbations but unstable against nonax- 
isymmetric perturbations (with m = 1 being the most un- 
stable mode). Differentially rotating PNS might be suscep- 
tible to that kink-type instability that in turn may limit their 
magnetic field amplitude. 

Criterion dHJ cannot be applied directly to fields under 
the influence of differential rotation. In a first step to under- 
stand the complicated interaction of toroidal magnetic fields 
and differential rotation we have studied a Taylor-Couette 
container with two corotating cylinders where the radial ro- 
tation law is hydrodynamically stable. An electric current 



1 see http://www.ioffe.rssi.ru/astro/conduct/ 

2 Other possible notations for are (3 = Rb : 



is flowing parallel to the rotation axis through the conduct- 
ing fluid, thus producing a nearly uniform toroidal magnetic 
field. It becomes unstable against nonaxisymmetric pertur- 
bations for non rotating cylinders but only for a rather strong 
magnetic field. If measured in terms of Hartmann numbers, 

Ha = J^, (9 ) 

with i? = \J (Rout — Rin)Rout, this is at about Ha = 150 
(see Fig. [TJ. In case of rotating cylinders without magnetic 
field the rotation law may be so flat that it is hydrodynam- 
ically stable. With magnetic fields it becomes always un- 
stable (Riidiger et al. 2007). One finds that differential ro- 
tation is strongly destabilizing for large magnetic Prandtl 
numbers which are characteristic for the PNS matter. One 
can thus expect that the toroidal magnetic fields induced by 
the differential rotation of PNS are limited by the described 
current-induced instability. 




50 100 150 
HARTMANN NUMBER 



200 



Fig. 1 The stability domain (hatched) for outer cylinder 
rotating with 50% of the inner cylinder (/Zfj = 0.5), the 
magnetic field is almost uniform (/is = 1), the perturba- 
tions are nonaxisymmetric (m = 1). Note the destabilizing 
action of high magnetic Prandtl numbers (here Pm = 100). 



In the present paper the stability problem for strong 
toroidal magnetic fields under the influence of differential 
rotation and the Hall effect is considered. We shall see that 
in this case even the sign of the toroidal field (with respect 
to the global rotation) will play an important role. 



3 The basic equations 

The basic state in the cylindrical system is Ur = U z = 
Br = B z = and 

B b 
B <j) = AR+-, U lj> = RO = aR+ -, (10) 

R R 
where a, b, A and B are constant values defined by 



a = Q\ 



A 



1 — i) 2 
Bin riifJ'B - ff) 



b — OinR in 



1 



1 — fj 



B — B in Ri, 



l-f) 2 ' 
1 - Pbt) 
I-77 2 ' 



(11) 
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Here is 

Rh 

V 



ft 



Ms 



Rout 

with i?; n and i? out as the radii, ft n and i? out the angular ve- 
locities, and B ln and B ou t as the azimuthal magnetic fields 
of the inner and the outer cylinders. The possible magnetic 
field solutions which do not decay are plotted in Fig. [2] 
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Fig. 2 The possible radial profiles of the toroidal mag- 
netic field between the two cylinders. The value of the in- 
tersection of each of the profiles with the right vertical axis 
is the corresponding ^-value. The profiles in the hatched 
domain are stable against axisymmetric perturbations while 
the cross-hatched area is also stable against nonaxisymmet- 
ric perturbations. The current-free solution B^ oc 1/R is 
given by the red line. 



We are interested now in the linear stability of the back- 
ground state ( [Tol l. In that case, the perturbed quantities of 
the system are given by 



u R: RQ + u<j>, 



B,i 



(13) 



As usual, the perturbations are developed in normal modes 
of the form 



F = F(R)exp(i(kz + rruj> + ut)). 



(14) 



Terms of the form (13[ and ([Pil l are inserted now into 
the the induction equation (O with the Hall effect included 
and linearized about the background state. The result is: 



°— - r?AB = E - pK 



with 

E R = 



R 



(imunB^ — imRQb R ) 



AB 't> , nu dhR 



-B, 



dun 
' dR 



BAku, 



Qb R + ikQRb z , 



(15) 



(16) 



(17) 



E, 



R 



(imUzB^ — \mRQb z 



(18) 



(12) and the Hall terms 



Hr 



1 

R 2 

-BAkRb 



dB 4 
dR 



■\kR bn + BchkmRb^ 



R 



B<f,m b z ,. 



(19) 



1 

R 2 



dfl rf 



-vmRb, 



dR 

-2B^\kRb^ — B^kmRbji 



db z 
B^imR—- 
dri 



BArnb- 



(20) 



H z = 



R 2 



,db 



R 



&{d& RbR+ dS" dR 

dBrh ^dbrh 
+2-J-Rb R + B^mR-^ 
dR dR 



dB d 



B A R 



dR 

db R 



iRbj 



dR 



+2BAmb ( t 



B^m 2 bi 



(21) 



The dimensionless numbers of the problem are the mag- 
netic Prandtl number (Pm), the Hartmann number (Ha) and 
the Reynolds number (Re), i.e. 

BinRo 



Pm 



Ha 



Re 



(22) 



where R = (i?; n (i? ut — -Rin)) 1 ^ 2 is the character- 
istic length scale, v the kinematic viscosity and rj the 
magnetic diffusivity. The magnetic Reynolds number is 
Rm = Pm Re and the Lundquist number is S = vPm Ha. 

We use Rq as a unit of length and R^ 1 as a unit of the 
wave number, rj/R as a unit of the perturbed velocity, Q m 
as a unit of angular velocity and ui, and B in as a unit of 
magnetic fields (basic and disturbed). 

In normalized quantities eq. ( Tl3T > may be cast in the fol- 
lowing form: 

iwRm b = D(b) + E — /3H (23) 
with 



E R — — (imBuji — i77ii?Rm Qb R 
R V 



RQ 



B'u R 

db R 
dR 4 



- \kBu z - 

dQ, 
R dR hR 



B- 



du 



R 



Rm 



dR 

Qb R + ikRQb 2 



E z = — \ Bu z - RmQRb; 

Here we have used the notations 
Q = Q- ln Q and B^ = B ln B. 
The diffusion terms are 

D R (b) 



(24) 

(25) 
(26) 

(27) 



d 2 b R 


m 2 


-k 2 b R 


dR 2 


1 db R 


2im 


b R 


R dR 




i? 2 ' 



(28) 
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d 2 b^ 
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171 u 


k%- 


dR 2 








1 d6 


2im 


b<f, 




R dR ' 


~ R 2 


and 








D z (b) 


d% z 
dR 2 


o 

m , 

— tOz — 

R 2 


k\- 



i? di?' 



(29) 



(30) 



In the same way the normalized momentum equation can be 
written as 



Re 



du 

di 



(UV)u+(uV)U = 

D(u) - VP + Ha 2 (curlB x b + curlb x B) , 
so that 



\uo Reu + Re G = D(u) 
with 

G R = mif2u R — 2Qu ( f > , 
G = (R 2 n)'^ +imOu 4 



VP + Ha L 



R 



G z = Qimu 
and 

im 



R 



Bbi 



4>< 



L, 



-{RB)'b R + ^Bb 4 



(31) 

(32) 

(33) 
(34) 
(35) 

(36) 

(37) 
(38) 



i 2 = +i—Bb z . 
H 

The perturbed flow as well as the perturbed magnetic field 
are source-free, i.e. 



dUR U R 771 

-dR + T +1 R^ +lkUz 







(39) 



.777, 



ikb z = 0. 



and 

db R bn_ 
dR R 

An appropriate set of ten boundary conditions is needed to 
solve the system. No-slip conditions as well as zero normal 
components for the velocity on the walls results in 

UR = u<j, = Uz = 0. (41) 

The boundary conditions for the magnetic field depend on 
the electrical properties of the walls. The tangential currents 
and the radial component of the magnetic field vanish on 
conducting walls hence 



dR + R =bR 



0. 



(42) 



These boundary conditions hold both for R = R m and for 

R = Rout- 



4 Results 

The equations have been solved for a simple model. The 
normalized gap width between the cylinders is 0.5 and the 
rotation law is rather flat approaching f2 oc R^ 1 hence 
Mfi = 0.5. The toroidal field in the gap is almost uni- 
form (fiB = 1) but it is not current-free. This field vi- 
olates ([8]l and is therefore Tayler-unstable with a critical 
Hartmann number of about 150 (Rudiger et al. 2007b). This 
instability is strongly modified by the differential rotation. 
The results are given in the Fig. [3] for the Hall parameters 
/3 = -0.01,0 and 0.01. The Hall parameter /3 and the 
magnetic Prandtl number are the free parameters of the sys- 
tem. Note, however, that due to (0 only \/Pm/3o is a phys- 
ical parameter in the definition of the Hall quantity $. As 
only the combination of /3oHa comes into the equations we 
can fix the sign of (3q and make the calculations for positive 
and negative Ha values or we can fix Ha as positive and use 
both signs of (3q. We prefer the second possibility so that the 
results for positive and negative /3q correspond to opposite 
magnetic field orientations. 

The solid line in Fig. [3] (bottom) is identical with the 
marginal limit between stability and instability in Fig. [T] 
We find that the system is destabilized by the rotation for 
high magnetic Prandtl numbers (Pm = 100). In contrast, 
for Pm = 1 the rotation stabilizes the flow (Fig. [3] top), 
which demonstrates very clearly the significant differences 
between the solutions with large and small magnetic Prandtl 
numbers. 

In all these cases, however, the Hall effect acts in the 
same direction. For positive (3 the stability domain is re- 
duced and for negative (3 the stability domain is increased. 
The stabilization (destabilization) of negative (positive) 
Hall (3 is a very common phenomenon of all the models. 
In other words, positive B^ (i.e. j3 > 0) lead to smaller crit- 
ical field amplitudes than negative B^, (i.e. j3 < 0). Hence, 
if indeed the nonaxisymmetric Tayler instability limits the 
strength of the induced toroidal fields B^, then the resulting 
amplitudes are different for different signs of B^ due to the 
action of the Hall effect. 



( 4 °) 4.1 Cell structure 



The cell structure of the neutrally stable modes is repre- 
sented by the resulting vertical wavenumber k. From the 
normalizations it follows the relation 



(43) 



Ro 



R'u 




for the vertical cell size in units of the gap width so that for 

fj = 0.5 

Sz 7T 

7? /T = T (44) 

Hence, for k ~ tt the cells are spherical while for k tt 
they are rather flat. Both possibilities are realized in the cal- 
culations. In Fig.|4]the results for Pm =100 are extended to 
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isymmetric Tayler vortices depends on that sign. If indeed 
realized in nature then the sign of the toroidal magnetic field 
(in relation to the rotation axis) can easily be read from the 
observations. 
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Fig. 4 The same as in Fig. [3] (bottom, Pm = 100 ). Sta- 
bility only exists at left from the curves. They are marked 
with the wavenumbers of the marginal instability. The num- 
bers show that also the cell structure strongly depends on 
the sign of the Hall effect. 
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Fig. 3 Tayler instability (to = 1) for various magnetic 
Prandtl numbers Pm and with Hall effect, [is = 1, /tin = 
0.5. The curves are labeled by the Hall parameter [3q. Note 
that positive (3 reduces the stability domain while negative 
(3 increases it. 



much higher values of the magnetic Reynolds number. The 
difference of the stability domains for different Hall param- 
eters grows even bigger as the rotation increases. The curves 
are marked with the corresponding wave number values. We 
find that for negative j3 the cells are spherical but they are 
rather flat for positive (3. Thus, not only the stability do- 
mains strongly differ for the Tayler instability for opposite 
signs of the Hall parameter but also the shape of the nonax- 



5 Growth rates 

As the Hall time is much longer than the rotation time the 
question arises whether the Hall effect strongly enhances 
the growth times of the Tayler instability. However, for the 
considered parameters the Hall effect hardly influences the 
growth rates of the Tayler instability (Fig. [3]). Neverthe- 
less, the Hall effect plays an important role for the stabil- 
ity map of the Tayler instability but not in the resulting lin- 
ear growth rates of the unstable disturbances. The growth 
rates are computed along vertical lines at Ha = 100 and 
Ha = 140 in Fig. [4] The curves are marked with the Hall 
parameter flo also including [3 = 0. The growth rates are 
given in units of the angular velocity of the inner cylinder, at 
the stability lines they vanish. A growth rate of 0.01 means 
an e-folding time of the instability of about 16 rotation pe- 
riods. For Pm = 100 this is the characteristic value for the 
Tayler instablity without Hall effect. This time is decreased 
by positive Hall effect and it is increased by negative Hall 
effect. For positive Hall effect the Tayler instability results 
as much faster than the Tayler instability for negative Hall 
effect. All the growth rates grow with growing Hartmann 
numbers. 

Even a weak Hall effect does not generally prolong the 
growth time of the Tayler instability which scales with the 
rotation time. In this case the Hall effect is only a modi- 
fication of another instability. Even if the Hall effect itself 
forms the instability (together with the differential rotation) 
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also then the ('shear-Hall') instability scales with the rota- 
tion rate and not with the rather long Hall time (Riidiger & 
Kitchatinov 2005). 

Another example for this phenomenon is given by the 
plane-wave solution of an a — ft dynamo. Both growth rate 
and cycle time of the most unstable mode of a linear oscil- 
lating ail dynamo with weak a-effect are mainly fixed by 
the basic rotation: 7/fi oc (w^/fi) 1 / 3 . Here 7 is the growth 
rate, u> v the dissipation frequency and il the basic rotation. 




4000 8000 
tnag. REYNOLDS NUMBER 




REYNOLDS NUMBER 



Fig. 5 The growth rates of the Tayler instability normal- 
ized with f2; n for m = 1 with and without Hall effect. The 
curves are marked by their Hall parameter (3 . The values 
belong to a vertical line for Ha = 100 (left) and Ha = 140 
(right) in Fig.gl Pm = 100, fi B = 1, (j,n = 0.5. 



6 Steeper magnetic profile 

The radial profile of the toroidal magnetic field used in Fig. 
|3] is rather smooth. Without detailed simulations one can- 
not not know the real profile. Hence, the computations rep- 
resented in Fig [3j are thus repeated for different magnetic 
field profiles for the most interesting case of high magnetic 
Prandtl number (Pm = 100). 

Figure [6] has been obtained for magnetic fields that in- 
crease outwards (/is = 3). According to the Tayler criterion 
dS) such profiles are highly destabilizing. Consequently, we 
find the critical Hartmann number for Re = one order 
of magnitude smaller than in Fig. [3] The opening of the two 
curves for (3 — ±0.01 is with about factor 2 for Rm ~ 1000 
very similar to the previous case. Again, generally the sta- 
bility domain for positive (3 is much smaller than for neg- 
ative (3. These basic findings do not depend on the actual 
Hartmann numbers for various magnetic profiles. Never- 
theless, we should underline that with the given parame- 
ters (p ~ 10 13 g/cm 3 ,i/ ~ 10 9 cm 2 /s,7y ~ 10 7 cm 2 /s) for 
Ha ~ 100 the maximal stable magnetic field is3xl0 12 G 
which value even grows for thinner layers. 

7 Asymmetry of the neutron star 
hemispheres 

Wardle (1999) has shown that due to the Hall effect the sta- 
bility properties of a differentially rotating MHD flow de- 
pend on the sign of the axial magnetic field. After our results 
the same is true for the azimuthal magnetic field. Moreover, 
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Fig. 6 Tayler instability (m = 1) for steep magnetic field 
(Pb = 3), for large magnetic Prandtl number (Pm = 100) 
and with Hall effect, = 0.5. The curves are labeled by 
their Hall parameter (3q. 



the critical magnetic field value above which the flow be- 
comes unstable can basically differ for different magnetic 
field orientations (Fig. [31 bottom). 

If the effect is strong enough this finding can have con- 
sequences. If in a PNS with differential rotation the toroidal 
field results from a poloidal field with dipolar symmetry 
then also the is antisymmetric with respect to the equa- 
tor. If the Tayler instability indeed determines the maximal 
field amplitudes then due to the Hall effect the amplitudes 
of the toroidal field in both hemispheres become differ- 
ent. Obviously, the Tayler-Hall instability produces an extra 
quadrupolar component of the originally produced toroidal 
fields with dipolar symmetry. It is thus unavoidable that the 
amplitudes of the induced toroidal field belts are different in 
both hemispheres. 

For Ha - 100, Pm - 100 and f3 Q ~ 10~ 2 taken from 
the bottom plot of Fig. fJJ we find for the Hall parameter 
{3 ~ 10 leading to ~ 10 13 G for the neutron star. This value 
is typical for pulsars so that the conclusion about different 
toroidal field values in both the hemispheres of a neutron 
star due to Tayler-Hall effects becomes realistic. 

On the other hand, strong magnetic fields suppress the 
heat transport in neutron stars (Schaaf 1988, 1990; Heyl & 
Hernquist 1998). The heat transport is blocked in the direc- 
tion perpendicular to the field lines so that the heat conduc- 
tivity tensor becomes anisotropic, i.e. 



Xij = Xx&ij +X2BiB 3 , 



(45) 



where \i represents the heat flux perpendicular to the field 
which is quenched by strong magnetic fields, hence (say) 
Xi oc 1/(1 + (3 2 ). With 



Xij 



(46) 



the heat flux remains finite along the field lines even for 

B -> 00. 
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The consequence of this magnetic-induced anisotropy 
of the heat flux tensor is a global inhomogeneity of the sur- 
face temperature (Geppert et al. 2006). If the latitudinal dis- 
tribution of the magnetic field is strictly symmetric or anti- 
symmetric with respect to the equator then the surface tem- 
perature results as equatorsymmetric. This is not true if for 
both hemispheres the magnetic amplitudes are differing (or, 
in other words, if the total magnetic field is a combination 
of a dipole and a quadrupole). Exactly this is the case if the 
toroidal magnetic field is produced by differential rotation 
under the presence of Tayler-Hall instability. If the differen- 
tial rotation of the neutron star disappears then the magnetic 
fields are frozen in so that the magnetic constellation is con- 
served for the time scales of the Ohmic decay (also modified 
by the Hall effect). We do thus expect the two half spheres 
of an isolated neutron star to be of different X-ray activity. 

Schwope et al. (2005) have indeed found an equatorial- 
asymmetric X-ray brightness analyzing XMM observations 
of the isolated neutron star RBS1223. The authors have as- 
sumed the existence of one bright "spot" in each of the 
hemispheres and found two temperature maxima of differ- 
ent strength (ratio e = 0.91). If this asymmetry effect is gen- 
eral for neutron stars then the interior magnetic fields must 
also be asymmetric with respect to the equator (dipole plus 
quadrupole) which can be explained with the Tayler-Hall 
scenario with differential rotation developed in the present 
paper. 
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